First Order Transitions and Multicritical Points in Weak Itinerant Ferromagnets 
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It is shown that the phase transition in low-T c clean itinerant ferromagnets is generically of first 
order, due to correlation effects that lead to a nonanalytic term in the free energy. A tricritical point 
separates the line of first order transitions from Heisenberg critical behavior at higher temperatures. 
Sufficiently strong quenched disorder suppresses the first order transition via the appearance of a 
critical endpoint. A semi-quantitative discussion is given in terms of recent experiments on MnSi, 
and predictions for other experiments are made. 
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The thermal paramagnet-to-ferromagnet transition at 
the Curie temperature Tc is usually regarded as a prime 
example of a second order phase transition. For materials 
with high Tc this is well established both experimentally 
and theoretically. Recently there has been a considerable 
interest in the corresponding quantum phase transition of 
itinerant electrons at zero temperature (T = 0), and in 
the related finite T properties of weak itinerant ferromag- 
nets, i.e. systems with a very low Tc- Experimentally, 
the transition in the weak ferromagnet MnSi has been 
tuned to different Tc by applying hydrostatic pressure 
PJ . Interestingly, the transition at low T was found to be 
of first order, while at higher transition temperatures it is 
of second order . The tricritical temperature that sep- 
arates the two types of transitions was found to roughly 
coincide with the location of a maximum in the magnetic 
susceptibility in the paramagnetic phase. Theoretically, 
it has been shown HQ that in a T — itinerant elec- 
tron system, soft modes that are unrelated to the criti- 
cal order parameter (OP) or magnetization fluctuations 
couple to the latter. This leads to an effective long-range 
interaction between the OP fluctuations. In disordered 
systems, the additional soft modes are the same 'diffu- 
sons' that cause the so-called weak-localization effects in 
paramagnetic metals p] . In clean systems there are anal- 
ogous, albeit weaker, effects that manifest themselves as 
corrections to Fermi liquid theory [|| . A Gaussian theory 
is sufficient to obtain the exact quantum critical behav- 
ior in the most interesting dimension, d — 3, for clean 
as well as for disordered systems (apart from logarithmic 
corrections in the clean case) pM]. 

In this Letter we show that at sufficiently low temper- 
atures the phase transition in itinerant ferromagnets is 
generically of first order. This surprising result is shown 
to be rooted in fundamental and universal many-body 
physics underlying the transition, viz. long- wavelength 
correlation effects, and hence to be independent of the 
band structure. This suggests that the behavior observed 
in MnSi is generic, and should also be present in other 



weak itinerant ferromagnets. We also make detailed pre- 
dictions about how quenched disorder suppresses the first 
order transition, which allows for decisive experimental 
checks of our theory. 

Let us start by deriving the functional form of the free 
energy of a bulk itinerant ferromagnet at finite T , and in 
the presence of quenched disorder that we parametrize by 
G = 1/cft, with cf the Fermi energy, and r the elastic 
mean-free time. The general Landau expansion of the 
free energy F as a function of the magnetic moment m 
in an approximation that neglects OP fluctuations is 



F = t m 2 + U4 m 4 + u e m 6 + 



la) 



The coefficients t, u^, Uq, etc. in this expansion can have 
nontrivial properties and contain important physics. A 
derivation from a microscopic theory shows that they are 
given as frequency-momentum integrals over correlation 
functions in a 'reference system' that depends on the na- 
ture of the underlying microscopic model Q . If the criti- 
cal magnetization fluctuations are the only soft modes in 
the system, then they are simply numbers. However, if in 
the process of deriving the Landau functional some other 
soft modes have been integrated out, then the coefficients 
will in general not exist, since they are represented as di- 
verging integrals over the soft modes. In Refs. || and jj] 
it was shown that in an itinerant electron system at T = 
there are indeed such soft modes. In the disordered case, 
these are the 'diffusons' mentioned above, with a disper- 
sion relation uj ~ fc 2 , and they lead to coefficients whose 
divergent parts have the form 
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Here A is a momentum cutoff, and all prefactors in the 
integrals have been omitted. In the clean case, the rele- 
vant soft modes are particle-hole excitations in the spin- 
triplet channel with a ballistic dispersion relation, w ~ fc. 
The resulting integrals are still divergent, although not 
as strongly as in the disordered case. It was shown in 
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Refs. j|, 0J that these divergent terms in the Landau 
expansion can be understood as an illegal expansion of a 
nonanalytic term in the free energy of the form 



/(to) = to / dkk 2 
Jo 



dw 
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In the disordered case, where x = 2, this follows explicitly 
from Eq. (3.6') of Ref. ||. In the clean case, an analo- 
gous treatment yields the same expression with x = 1. 
Notice the different sign of the dirty case compared to 
the clean one, which we will come back to below. Equa- 
tion yields /(to) cx to 5 / 2 and /(to) oc to 4 In to in the 
disordered and clean cases, respectively. In either case 
the resulting singularity is protected by the magnetiza- 
tion, which gives the soft modes a mass. The leading 
effect of T 7^ is adequately represented by replacing 
lj — > u> + T . In addition, in the presence of disorder the 
ballistic modes in the clean case obtain a mass propor- 
tional to 1/r, so the appropriate generalization of Eq. 
(|J) for the clean case (x = 1) to finite temperature and 
disorder is obtained by the replacement to — > uj + T+1/t. 
Doing the integrals, and adding the usual terms of order 
to 2 and to 4 , we obtain a free energy of the form 

F = t to 2 + G {N F T t ) to 4 [m 2 + (<xT) 2 ] ~ 3/4 

+v to 4 In (to 2 + (T + (3Gf) +um i + 0(m 6 ) , (3) 

where T t is an effective spin-triplet interaction amplitude 
H made dimensionless by means of a density of states at 
the Fermi level, N-p. If we measure F, to, and T in terms 
of a microscopic energy, e.g. ep, then t, v, and u are all 
dimensionless. v is quadratic in T t 0J. t = 1—NpU is the 
dimensionless distance from the critical point. It depends 
on the physical spin-triplet interaction amplitude U, with 
N-p U w 1 in a ferromagnetic or nearly ferromagnetic sys- 
tem, while r t above is an effective interaction amplitude 
with NpT t < 1. r t is expected to be relatively larger in 
strongly correlated systems. Finally, a and (3 are param- 
eters that measure the relative strengths of the tempera- 
ture and the disorder dependence, respectively, in the two 
nonanalytic terms. They are numbers of order unity, and 
like u and v they are non-universal. Equation (|^) pro- 
vides a functional form of the free energy that correctly 
describes the leading nonanalytic m-dependence for both 
clean and disordered systems, as well as the leading tem- 
perature cutoff for either term and the leading disorder 
cutoff for the clean nonanalyticity. 

The sign of v merits some attention. Perturbation the- 
ory to second order in T t yields v > BM . Further, v > 
indicates a decrease of the effective Stoner coupling con- 
stant I due to correlation effects: / is a homogeneous 
spin susceptibility, v > means that this susceptibility 
increases as the wavenumber increases from zero jq] , and 
correlation effects decrease with increasing wavenumber. 
It is well known that correlation effects in general de- 
crease / [01, and v > is consistent with that. Ref. Hi 
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FIG. 1. Phase diagram at T = for u = 1, v = 0.5, 
N-pTt = 0.5, a — /3 — 1, showing a second order transition 
(dashed line), and a first order transition (solid line). 



has given some possible mechanisms for v to be negative 
at least in some materials, and showed that in this case 
the ferromagnetic transition is always of second order. 
However, the generic case is v > 0, which we will now 
discuss. 

We first consider the case T = Q. The transition in 
the clean system, G = 0, is then of first order, since 
to 4 In to < for small m. Upon disordering the system, 
G > 0, the negative term is no longer the leading one at 
t = 0. For small values of G, the transition remains first 
order. However, for G exceeding a value G ce the first or- 
der transition occurs only at t < 0, and it is pre-empted 
by a second order transition. Since the negative term is 
only the third term in an m-expansion of F, the multi- 
critical point where the nature of the transition changes 
is a critical endpoint (CEP) p0| . The phase diagram in 
the G-t plane is shown in Fig7|l[. For G ce < G < G c , the 
second order transition at t = is followed by a second 
transition, the second one being of first order, to a state 
with a larger magnetization. The line of first order tran- 
sitions ends in a critical point (CP) at a disorder value 
G c , where the two minima in the free energy merge. 

Before we consider T > 0, let us discuss this re- 
sult and the validity of our conclusions. To facilitate 
an analytic discussion, we put = 0. We then have 
F = tin 2 + G(N F T t )m 5 / 2 + 2vm i \nm + urn 4 . At G = 
there is a first order transition at t = uexp[— (1 + u/v)], 
and the magnetization at the transition has a value 
to = cxp[— (1 + u/v)/2}. Notice that the nonanalytic 
term is the leading one in F after the tm 2 term, and 
that we know the functional form of F exactly up to 
0(m ). As long as u/v >> 1, m is exponentially 
small at the transition. For small v, our Landau ex- 
pansion is therefore controlled in the sense that terms 
of 0(to 6 ) and higher would have to have exponentially 
large coefficients in order to change our results. For 
G > G ce = (4u/3JV F r t )exp[-(l + 3u/lv)], the first 
order transition is pre-empted by a second order one. 
At the CEP, the magnetic moment has the value to = 



exp[-(2/3 + u/2v)} 



-1/6 



m(G — 0). Allowing for 



f3 = 0(1) 0, and repeating the calculation numerically, 
leads only to minor quantitative changes of these results. 
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FIG. 2. Phase diagrams for u = /3 = 1, v = a = NpTt = 0.5 showing first (solid) and second (dashed) order transitions. 



At T > 0, the free energy is an analytic function of 
m, but for small T the coefficients in an m-expansion be- 
come very large. Our remarks about the validity of our 
truncated Landau expansion therefore still apply, i.e., at 
< k^T « e-p, our theory contains the most important 
terms to every order in an expansion in powers of m 2 . 
Let us first consider the clean system, G = 0. There is a 
tricritical point (TCP) at T tc = exp(—u/2v), with a first 
order transition for T < T tc , and a line of Heisenberg crit- 
ical points for T > T tc . To describe the (conventional) 
tricritical behavior in d = 3 our mean-field theory is suf- 
ficient (apart from logarithmic corrections) PH , for the 
critical behavior at T > T tc it is of course not. 

For the suppression of the first order transition by dis- 
order at T > we find two different possibilities, de- 
pending on the value of the parameter a. For small a 
(a < 1.5 with our choice of the remaining parameters, 
Fig. |), the TCP is replaced by a CEP for G larger than 
some G tce < G ce . At G = G ce , the CEP reaches T = 0, 
and for larger values of G the transition is of second order 
for all T. At small T, it is followed by a first order tran- 
sition. The line of first order transitions ends in a critical 
point, and disappears only for G — G c . For larger values 
of a (Fig. ||), the TCP persists for a range of disorder 
larger than G ce . The first order transition first gets pre- 
empted in a temperature window between two CEPs. At 
G = G ce , the lower CEP reaches T = 0, while the TCP 
at higher temperature survives. With further increas- 
ing disorder, two CPs appear in the ordered phase, and 
the remaining CEP gets replaced by a TCP. Finally, the 
two TCPs merge, and the remaining CP reaches T = 0, 
eliminating the last temperature regions with first order 
transitions. Notice that the interesting features of these 
phase diagrams do not depend on the logarithm in Eq. 

; similar features are obtained in standard phenomeno- 
logical Landau expansions with a negative coefficient of 
the third term p2| . We stress again, however, that in our 
case the expansion is controlled, and that we have a defi- 
nite physical mechanism for the appearance of a negative 
term, in contrast to purely phenomenological theories. 

We now turn to a discussion of the available exper- 
imental information on this subject. MnSi has a low 
Tc (~ 30 K) under ambient pressure, and T c can be 



driven to zero by a hydrostatic pressure p c ~ 15kbar. 
k^T ' I 'e-p << I always, and T is low enough to suppress 
phase breaking processes, so the quantum critical behav- 
ior is easily accessible experimentally. This system has 
been studied in detail by Pfleiderer et al. jjj These au- 
thors found from susceptibility measurements that the 
transition turns first order at a T c of about 12 K. The line 
of second order transitions was found to scale with pres- 
sure like T c cx (p c — p) 3 / 4 , while in the first order regime 
the transition temperature varies like T± cx (p — pc) 1 / 2 . 
The scaling of T c with pressure was explained by a scal- 
ing analysis based on the self-consistently renormalized 
(SCR) theory of Moriya and Kawabata jl3| , assuming a 
dynamical exponent z = 3. The first order transition at 
low T was attributed in Ref. Q to a sharp structure in 
the density of states at the Fermi level. 

Let us look at the experiment in the light of the above 
discussion. In Ref. Q it was shown that the quantum 
phase transition in d = 3 is indeed correctly described 
by SCR theory, apart from logarithmic corrections that 
would be very difficult to detect experimentally, and that 
the dynamical critical exponent in d = 3 is z = 3. The 
analysis of Ref. Q was therefore adequate, and in par- 
ticular the quantum-to-classical crossover exponent <f>, 
which determines the behavior of the critical temperature 
as a function of t through the relation T c oc tP, has a value 
<f> = 3/4. If one makes the plausible assumption that t 
depends linearly on the hydrostatic pressure, at least for 
small t, then this is in agreement with both the exper- 
imental finding and the analysis in Ref. Q. As for the 
pressure dependence of T\ , one of the temperature scales 
in the problem is the Fermi liquid temperature scale [Q, 
which arises from a quadratic T-dependence of t. Since 
the first order transition is determined by the condition 
t(T±) — const., we immediately get T\ cx y/p c — p, where 
we again assume a linear relation between p and t. 

We finally discuss the observation Q that the tricritical 
temperature roughly coincides with a minimum of the 
inverse magnetic susceptibility x _1 m the paramagnetic 
phase. In d dimensions, the leading T-dependence of the 
paramagnetic susceptibility is of the form || 

X /2N F = l + 2i d T 2 T d ~ 3 ~u d T 2 , (4) 
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FIG. 3. Same as Fig. § but for a = 2. 



In d = 3, the nonanalyticity is of the form T 2 \nT. A 
calculation of v$ to second order in T t revealed || that 
to that order, V3 = 0, in agreement with prior results 
from Fermi liquid theory |l5) . Ref. || also discussed that 
there are reasons to believe that the exact value of V3 = 
may be nonzero. If we assume that this is the case, then 
we obtain a minimum in at a temperature T nim — 
cxp(— U3/2V3 — 1/2). Since the nonanalyticities in F and 
X are manifestations of the same singularity, one expects 
U3 w u and V3 m v, so that T mlu w T tc . While this 
provides a possible explanation for the observation, we 
stress the speculative nature of the above considerations 
due to the theoretical uncertainty about a nonanalytic T 
dependence of x m d = 3. 

Our theory thus provides us with a complete explana- 
tion for the nature of the transitions observed in MnSi, 
and in particular for the existence of a first order tran- 
sition at low T, which in Ref. jl) was attributed to a 
band structure feature characteristic of MnSi. While 
this feature may well be sufficient to make the transi- 
tion in MnSi of first order, the present theory leads to 
the surprising prediction that the first order transition is 
generic, and thus should be present in other weak itiner- 
ant ferromagnets as well. Our theory further predicts in 
detail how the first order transition will be suppressed by 
quenched disorder. Observations of such a suppression, 
or lack thereof, would be very interesting for corrobo- 
rating or refuting the theory. Semi-quantitativcly, the 
theory predicts that the T region that shows a first or- 
der transition will be largest for strongly correlated sys- 
tems. Conversely, since the dependence of the tricritical 
temperature on the system parameters is exponential, in 
some, or even many, systems the first order transition 
may take place only at very low temperatures. This may 
explain why in ZrZn2 no first order transition has been 
observed Q|, although the experiment does not seem to 
rule out a weakly first order transition [Q . 
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